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THE EFFECT OF AN ISOTROPIC NON-EQUILIBRIUM 
PLASMA ON ELECTRON TEMPERATURE MEASUREMENTS 

Robert F. Benson 
Walter R. Hoegy 

Laboratory for Planetary Atmospheres 
ABSTRACT 

The electron temperatures which would be determined (vising 
the conventional single-temperature analysis) by the electro- 
static pi’obe, the diffuse resonance, and the radar backscatter 
techniques in an isotropic two temperature plasma are pre- 
sented. Plasma models corresponding to the addition of a 
minor component of energetic electrons, and models cor- 
responding to a process that cools a fraction of the ionospheric 
electrons are considered. The diffuse resonance temperature 
is found to lie between the probe and radar backscatter tempera- 
tures. The isotropic models corresponding to the addition of 
energetic electrons cannot support the reported discrepancies 
between radio wave and probe electron temperature measure- 
ments. Temperature differences similar to the observed 
differences can be produced by models with a fraction of the 
electrons at a temperature cooler than that of the main com- 
ponent of electrons. These models, however, are difficult to 
explain in terms of our present understanding of the ionospheric 
plasma. 
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INTRODUCTION 


A discrepancy has been observed between the electron temperature de- 
duced from in situ electrostatic probe measurements and from ground-based 
radar backscatter measurements; the probe temperature T p being greater than 
the backscatter temperature T b (Hanson et al. , 1969; Carlson and Sayers, 1970). 
A similar discrepancy has been observed in a comparison of simultaneous in 
situ Alouette 2 electron temperature measurements based on the rf sounder 
stimulated plasma resonances and the cylindrical probe experiment (Oya and 
Benson, 1972). The resonance measurements are based on the diffuse reso- 
nance which occurs at the frequency f DJ , between f H and 2f H (Oya, 1970), where 
f H is the electron cyclotron frequency; the probe measurements are based on 
the experiment of Findlay and Brace (1969). Again, the radio wave temperature, 
i. e. , the diffuse resonance electron temperature T f , is less than the probe 
temperature T p . 

An assumption inherent in the determination of T p , T b , and T| is that 
the electron velocity distribution function f(v) is Maxwellian. The purpose of 
the present paper is to determine the effect of an isotropic non- Maxwellian 
velocity distribution on the three temperature measuring techniques. For 
illustrative purposes we represent the distorted velocity distribution by a 
simple, isotropic two temperature distribution - a superposition of two 
Maxwellian distributions. Hoegy (1971) applied this distribution to a study of 
the probe-radar temperature difference for distortions in f(v) at high energies 
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(photoelectron distortion) and at low energies (thermal distortion). It was con- 
cluded that significant temperature differences could be produced by a thermal ■ 
distortion but not by a distortion due to realistic photoelectron populations. 

This paper represents an extension of the work of Hoegy (1971) to include the •' 
diffuse resonance technique. 

Basically, we pose the following question: if an isotropic two temperature 
plasma is present, but the data are interpreted in terms of a one temperature 
Maxwellian plasma, what temperature will be computed? A description of the 
method for obtaining T p and T b for such a plasma has been given by Hoegy 
(1971); the procedure for obtaining T f is given in the Appendix. 
ILLUSTRATIVE EXAMPLES 

The effect of an isotropic distorted f(v) on the temperatures deduced by 
the probe, radar, and diffuse resonance techniques is illustrated in Figure 1 
for a wide range of conditions. In all cases the bulk of the electrons are con- 
sidered to be at 3,000°K and the perturbing electron component extends over a 
range of temperatures above (top row) and below (bottom row) this value. By 
convention, T, denotes the temperature of the cooler electron component re- 
gardless of its role as the major or minor constituent. The three columns of 
figures from left to right in Figure 1 represent perturbing electron populations 
of 10%, 20%, and 30%, respectively. 
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The following trends are evident from Figure 1: .. . . 

(1) T, is between T p and all .cases^, , , . , 

(2) A larger perturbation in the j^robe temperature ca.n be obtained when 
a fraction of the electrons are warmer than the main component 
rather than cooler than the main component; the converse is true 

in the case of the diffuse resonance and radar backscatter tempera- 
tures (compare Figure la with Id, lb with le., or lc with If). 

(3) All three techniques give very similar results when warmer electrons 
(at T 2 ) are present up to T 2 * 10, 000°K; when electrons with 

T a > 10,000°K are present the radio wave temperatures ti (Tf and 
T b ) are greater than the probe temperature (see Figures la, lb, 
and 2c). ; ... 

(4) The radio wave temperatures are always less than the probe tem- 
peratures (T b < T, oj < T p ) when cooler electrons are present (see 
Figures Id, le, and If). 

The following comments should be kept in mind when interpreting Figure 1. 
(1) The actual temperature values for T p and T f[)] depend on the simulated ex- 
perimental conditions - the voltage interval for the probe and the ambient plasma 
conditions for the diffuse resonance. (2) There is a double-valued behavior of 
T p as a function of the perturbing temperature because the perturbing current 
can be recognized and rejected in the temperature reduction process when the 
perturbing temperature differs by a large amount from the temperature of the 
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major component of electrons. (3) The peak observed in T, 0) in the top row of 
Figure 1 corresponds to the conditions giving the maximum distortion in the two 
temperature dispersion curves (see Figure A1 of the Appendix). 

DISCUSSION 

Hanson et al. (1969) observed T p to be as much as 70% greater than t., 
and Oya and Benson (1971) found T p to be typically 20% to 30% greater than 
T, . It is clear that these discrepancies cannot be explained in terms of an 
isotropic plasma with a component of hot electrons present (see comment 3 of 
the previous section). 

If an isotropic plasma exists with a fraction of the electrons at a tempera- 
ture cooler than the main component, then it is possible to produce a temperature 
difference resembling the observed discrepancies (see comment 4 of the previous 
section). For example, consider the case of a 10% perturbation (Figure Id) with 
T, = 1, 200°K; the two temperature model then gives the following values: 

T p = 2930°K, Tf Dj - 2740°K, and T b = 2610°K. The calculated ratios for 
T p/ T fDi anc * T p/ T b (1*07 and 1.12, respectively) are smaller than the observed 
ratios, but larger values can be obtained by using a lower T, and/or a higher 
minor constituent population (see Figures Id, le, and If). 

In order to maintain such large distortions of the low energy electron 
population, there must exist heating and cooling processes with rates comparable 
in magnitude to the electron-electron energy exchange rate. Since this time is 
three orders of magnitude smaller than the electron- ion or electron-neutral 
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energy exchange time, the main effect of Coulomb collisions with cooler particles 
is to modify the average energy of the electrons rather than to modify their 
Maxwellian velocity distribution. Other processes, however, may be important 
enough to modify this picture-such as non-local processes or wave-particle 
interactions which could enhance the selective heating and cooling mechanisms 
at the expense of the electron-electron equilibration. This possibility merits 
future investigation. 

The present work has treated the isotropic case. Recent observations by 
Clark et al. (1973) indicate that significant anisotropies, with the electron tern- 

— V 

perature parallel to the earth's magnetic field vector B greater than the tem- 
perature perpendicular to B, are common in the upper ionosphere. Such 
anisotropies could contribute to the observed discrepancy between radio wave 
and probe measurements (for example, T <0j corresponds to electron motion 
perpendicular to B, whereas T p is not very sensitive to direction). Thus, in 
addition to considering mechanisms which could produce a distortion in the 
isotropic distribution of electron velocities, future work should consider the 
anisotropic case. 

APPENDIX 

t, d1 IN A TWO TEMPERATURE PLASMA 
In order to determine T| qi in a two temperature plasma it is necessary 
to derive the two temperature dispersion relation for plasma waves propagating 
perpendicular to the earth's magnetic field vector B. (The diffuse resonance 


5 



experiment is interpreted in terms of perpendicular propagation (Oya and 
Benson, 1972).) 

The dispersion relation for longitudinal plasma waves with the wave vector 
k perpendicular to B is obtained by setting the expression for the dielectric 
constant ( e long .) A e( ! ua l to zero, ie. , 


= ( €| ) =l + ^ei 


+ 00 


2ttv , dv 
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- oo 


is the distribution function for electron motion perpendicular to B, e and m are 
the electron charge and mass, v„ and v ± are the components of electron velocity 
v parallel and perpendicular to B, co=27rf, J2=27rf H , k= lk| , f (v) is the electron 
velocity distribution function, and J„ is the Bessel function of order n (e.g. , 

see the review by Dougherty and Watson (1971)). Inserting the isotropic two- 
temperature velocity distribution function 



(A2) 


into Al (where p represents the population of electrons at T a ) gives the two 
temperature dispersion relation 
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The standard single-temperature dispersion relation is obtained by setting 
p=0 or 1 or by setting T, =T, in (A3). Note: When very energetic electrons 
are considered, large X 7 values are encountered and the second summation 
term in (A3) converges very slowly. A much faster convergence is obtained by 
making use of the identity 


+ OO 

Y J e- x I„ (X) = 1 


in order to make the following replacement 



The dispersion curves based on (A3) are presented in Figure A1 in the 
form of the frequency f, normalized by the electron cyclotron frequency f H , 
versus the wave number k normalized by 1/R, where R represents the electron 

/ Xl/^ 

cyclotron radius. In Figures la-lc, R is expressed as R=R(T,) = lKT,/m(27rf H ) ] 
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where T, is the temperature of the major component; in Figures ld-lf it is ex- 
pressed as R=R(T 2 ) = ^KT 2 /m(27rf H ) 2 y/ 2 where T 2 is the temperature of the 
major component. The first case (Figures l(a)-l(c)) corresponds to the ad- 
dition of energetic photoelectrons, while the second case (Figures l(d)-l(f)) 
corresponds to a process that cools the ionospheric electrons. The solid line 
in each figure corresponds to the single temperature dispersion relation that 
was used by Oya and Benson (1972) in the determination of T, . 

T D1 

The procedure for determining T, 01 from a single-temperature dispersion 
relation in a two temperature plasma is as follows: Select the two temperature 
parameters p, T,, and T 2 , and the ionospheric parameters cu N , fi, k, that 
correspond to the conditions of the observations. Then determine the nor- 
malized diffuse resonance frequency from the two temperature dispersion re- 
lation (A3). Finally, use the standard single temperature dispersion relation 
to obtain A= (kR) 2 and hence T fpi = Xmfi 2 /Kk 2 . 
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FIGURE CAPTIONS 


Figure 1. The electron temperature T p , T b , T foi deduced from a single 
temperature analysis using a two temperature plasma for the 
electrostatic probe, the radar backscatter, and the diffuse reso- 
nance techniques, respectively. T versus T 2 with T 2 >T, = 3,000°K 
in the top row, and T versus T, with T,<T 2 = 3,000°K in the 
bottom row. The left column corresponds to the condition when 
10% of the electrons deviate from the model temperature of 3,000°K, 
the center column corresponds to a 20% perturbation, and the right 
column to a 30% perturbation. Note the abscissa scale change be- 
tween the top and bottom rows. In each column, the bottom figure 
represents an enlarged extension of the top figure. to the temperature 
range below 3,000°K. The T p values correspond to the 0.3 to 1.3v 
interval of the probe current-voltage curve; the T, 0( values cor- 
respond to the plasma conditions given in the caption of Figure Al. 

Figure Al. Dispersion curves based on the two-temperature dispersion equation 
(A3) for plasma models with T 2 >T, = 3,000°K (a), (b), and (c); and 
for models with T, < T 2 = 3,000°K (d), (e), and (f). The parameter 
p represents the fraction of electrons at T 2 . The solid curves 
(p=0 and 1) correspond to the single temperature dispersion curve. 
The plasma conditions used in these models correspond to a typical 
observation from the Ottawa data sample of Oya and Benson (1972) 
2 tt - 3 . 475 MHz and f H =T.178 MHz so that f N /f H =2. 95; the 
corresponding value for kR when the temperature is taken as 
.3, 000°K is 0.97). 
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